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ABSTRACT 

We outline our techniques to characterise photospheric granulation as an astrophysical noise source. 
A four component parameterisation of granulation is developed that can be used to reconstruct stel- 
lar line asymmetries and radial velocity shifts due to photospheric convective motions. The four 
components are made up of absorption line profiles calculated for granules, magnetic intergranular 
lanes, non-magnetic intergranular lanes, and magnetic bright points at disc centre. These compo- 
nents are constructed by averaging Fe I 6302 A magnetically sensitive absorption line profiles output 
from detailed radiative transport calculations of the solar photosphere. Each of the four categories 
adopted are based on magnetic field and continuum intensity limits determined from examining three- 
dimensional magnetohydrodynamic simulations with an average magnetic flux of 200 G. Using these 
four component line profiles we accurately reconstruct granulation profiles, produced from modelling 
12 x 12 Mm 2 areas on the solar surface, to within ~ ± 20 cm s _1 on a ~ 100 m s" 1 granulation 
signal. We have also successfully reconstructed granulation profiles from a 50 G simulation using the 
parameterised line profiles from the 200 G average magnetic field simulation. This test demonstrates 
applicability of the characterisation to a range of magnetic stellar activity levels. 
Subject headings: Line: profiles - Planets and satellites: detection - Sun: granulation - Stars: activity 
- Stars: low-mass - Techniques: radial velocities 



1. INTRODUCTION 

One of the consequences of plasma motions in the outer 
layers of low-mass stars with convective envelopes is ra- 
dial velocity (RV) shifts due to variable stellar line profile 
asymmetries, known as astrophysical noise (or stellar 'jit- 
ter'). This can pose a major problem for planet hunters 
as RV follow-up is mandatory for most planet confirma- 
tion and characterisation. As these net RV shifts, pro- 
duced by photospheric convective m otions, are on the m 
s _1 level ()Schriiver fe Zwaanl I2000D . stellar granulation 
will severely impact the confirmation of low-mass plan- 
ets around solar-like stars. Additionally, the promise of 
cm s _1 RV precision from spectrographs like Espresso 
and Codex makes it ever more pertinent that we under- 
stand granulation both in terms of the physical processes 
involved as well as an important source of astrophysical 
n oise. 

IMeunier et al.l ()2010D conclude that the limiting factor 
for detecting Earth-like planets in the habitable zone is, 
in fact, the convective contribution to the RVs and not 
the larger amplitude contribution of spots and plages. 
IMeunier et al.l (|2010f) go on to affirm that the amplitude 
of convective noise in solar-like stars, must decrease by 
one order of magnitude or more to make possible the 
detection of low-mass planets. 

To date, there have been only a handful of attempts 
to decrease the effects of granulation for the benefit of 
exopla net detection and confirmation. iDumusaue et all 
(|20ll propose adjustments to the number of measure- 
ments and exposure time in order to average out some 
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of the effects of granulation. Alth ough this method has 
had success (see IPepe et al.l [2011f ). it is both time and 
cost intensive, and do es not provide inform ation on the 
nature of granulation. lAigrain et al.l (|2012h take a differ- 
ent approach and use photometry to estimate RV vari- 
ations due to stellar activity, including a parameter to 
estimate the perturbation due to convection. They suc- 
cessfully demonstrate that they can reproduce variations 
well below the m s _1 level. However, as low-mass planets 
around solar-like stars require cm s _1 precision, there is 
still an urgent need to further characterise and remove 
the RV variations of surface convection. 

Although there have been few attempts to remove or 
reduce the effects of granulation, there have been numer- 
ous observations and theories devoted to understanding 
the nature of granulation (for a compr ehensive review 
see iNordlund et al.l (|2009D . ISteinl (|2012D . and references 
therein). We we re particularly inspired by the four com- 
ponent model of IDravinsI (|1990f) . used to reconstruct the 
observed asymmetri es in stel l ar line profiles due to gran- 
ulation. Similar to IDravinsI (|1990D . our aim is to mea- 
sure and understand the RV signal produced by convec- 
tive motions integrated across a stellar disc. However, 
our main goal is to simulate numerous granulation obser- 
vations and ultimately develop a statistically significant 
noise removal method. 

In this paper, the first in a series, we outline a tech- 
nique to characterise granulation as an astrophysical 
noise source. The backbone of this characterisation 
is a state-of-the-art three-dimensional magnetohydrody- 
namic (MHD) solar simulation, coupled with detailed 
wavelength-dependent radiative transfer. Due to the 
time-intensive nature of detailed radiative diagnostics, 
producing enough realistic granulation patterns to cover 
an entire stellar disc with this method is not feasible. 
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As such, we use the output from a short time-series of 
simulations to parameterise the granulation signal. 

In Section [2] we outline our procedure to separate 
absorption line profiles from solar simulations into four 
categories based on continuum intensity and magnetic 
field strength to categorise the different physical com- 
ponents present in photospheric magneto-convection. In 
Section [3l we use these components to reconstruct the 
absorption line profiles from the simulated models. In 
Section |4] we show that this parameterisation can be suc- 
cessfully applied to other magnetic field strengths. We 
discuss our procedure to generate new realistic granula- 
tion line profiles in Section [5] and present the success of 
our parameterisation as well as its current limitations in 
Section [6] 

2. CHARACTERISATION OF CONVECTION 
COMPONENTS 

The radiative MHD simulations of solar magneto- 
convection reproduce the properties of granulation ob- 
served at the solar surface in detail and for a wide range 
of magnetic activity levels (for a comparison of simula- 
tion and observation see Fig ure D] and for furthe r com- 
parisons with observables see Shcl vag et al.l[200l . How- 
ever, there is no established to method to remove the 
granulation signal from observational data. We use such 
simulations, examining the Sun as a star in this paper, 
to develop a parameterised model of granulation that 
will eventually be used to remove convective noise from 
observations and can later be applied to other spectral 
types. Applicability to other spectral types is essential 
to developing a robust noise removal technique for exo- 
planet surveys due to the variety of planet-hosting stars. 
The simulations also provide many advantages towards 
developing such a technique including controllability of 
magnetic flux, seeing/resolution, and rotation. 

We simul ate solar surface gr anulation using the MU- 
RaM code (jVogler et al.l 120051 ) . The code solves large- 
eddy radiative 3D MHD equations on a Cartesian grid, 
and employs a fourth-order central difference scheme to 
calculate spatial derivatives. The numerical domain has 
a physical size of 12 x 12 Mm 2 in the horizontal direc- 
tion and 1.4 Mm in the vertical d irection and is resolv ed 
by 480 x 480 x 100 grid cells fsee lShelvag eFaIl[20lTI for 
details) . 

We first introduce a uniform vertical 200 G magnetic 
field and produce a sequence of 190 snapshots of solar 
magnetic granulation, with a ~ 30 second cadence. The 
sequence covers approximately 80 minutes of physical 
time, corresponding to ~ 10 - 20 granular lifetimes. This 
does not include the first 20 minutes of the simulation, 
which is discarded, in order to allow the large amplitude 
oscillation, induced from the introduction of magnetic 
field, to dampen and the m odel to stabilise (| Vogler et al.l 
l200alVogler fc Schussle"rll2007l) . During this first part of 
the simulation, the magnetic field is advected into the in- 
tergranular lanes by the convective plasma motions, and 
the maximum field strength rises from its initial value of 
200 G to a few kilogauss in the intergranular lanes. The 
200 G field is sufficiently strong to ensure each physical 
aspect of convection found in the quiet Sun is present in 
the simulations. 

We also produce simulations in a similar manner for 
a uniform vertical 50 G field. The cadence again is ~ 



30 seconds, and this simulation covers approximately 45 
minutes of physical time. A similar oscillation in the 
50 G model is present throughout the simulation since 
the lower magnetic field does not provide as much damp- 
ing compared to the 200 G model. This weak magnetic 
field simulation provides a valuable test case for our pa- 
rameterised model. An accurate reconstruction of the 
50 G model from the 200 G components demonstrates 
the robustness of our parameterisation across magnetic 
field strengths. 

The calculated models contain basic solar photospheric 
plasma parameters such as density, pressure, tempera- 
ture, magnetic field strengths, and velocity. We perform 
the radiative transport diagnostics with these parame- 
ters, using a modified version o f the STOPRO code (as 
described in lShelvag et al.l l2007). for the simulated pho- 
tosphere with the Fe I 6302 A magnetically sensitive ab- 
sorption line in local thermodynamic equilibrium (LTE) 
approximation. This line is widely used in both solar ob- 
servations and simulations as a diagnostic tool fo r mag- 
netic field and temperature ([Shelvag et al.l [200 7). The 
full Stokes vector is calculated for each pixel in the snap- 
shots throughout the simulation; for this study we use 
only the / component of Stokes vector. The wavelength 
range used for calculation is ± 0.3 A and is resolved by 
400 wavelength points. 

2.1. Four Component Model 

The premise for our parameterisation of granulation 
line profiles is that they are a summation of absorption 
line profiles corresponding to each physical component of 
solar surface granulation. As such, we began our param- 
eterisation by examining the different physical compo- 
nents of photospheric magneto-convection in the 200 G 
simulations. Granules are non-magnetic, hot, bright bub- 
bles of rising gas. The plasma in granules cools, darkens 
and falls back down into the intergranular lanes. The 
magnetic bright points (MBPs) present within the dark 
intergranular lanes are due to the enhanced continuum 
intensity and decreased radiation absorption in evacu- 
ated, highly magnetic flux tubes. Thus, a four com- 
ponent parameterisation, which divides the granulation 
pattern into granules, magnetic intergranular lanes, non- 
magnetic intergranular lanes and MBPs is a natural step 
to move forward. 

The simulation provides information on magnetic field 
and continuum intensity (an example of which is shown 
in the left panel of Figure [1]) , alongside the Stokes pro- 
files. Since the four components differ inherently in 
brightness and magnetic field, it is logical to impose cut- 
off limits based on continuum intensity and magnetic 
field (measured at the continuum formation height) to 
separate each component within the simulation. As a 
result, for each snapshot we average together the Stokes 
I profiles that correspond to each category to create our 
component profiles. Thus, for each of the 190 snapshots 
in the time sequence there are four line profiles, each rep- 
resenting a component within the granulation pattern. 

To determine the optimal magnetic field and contin- 
uum intensity limits for the four components we explored 
the parameter space from 100 - 1200 G and from 0.8 
to 1.2 times the average continuum intensity. The final 
selection of limits, 0.9 average continuum intensity and 
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Fig. 1. — An example of our simulation of solar granulation (left) in the white-light continuum (417 nm) compared to observations (right). 
The observations were taken with the Rapid Oscillation in the Solar Atmosphere (ROSA) imaging system on the Dunn Solar Telescope. 
The contrast in the observations is reduced due to seeing effects. 



1000 G, was based on how well the four profile compo- 
nents could reconstruct the original simulated line pro- 
files; the procedure of reconstruction is discussed in Sec- 
tion [3] The four profiles used to characterise the differ- 
ent components of photospheric magneto-convection are 
shown in Figure [5J The breadth about each component's 
respective means is representative of the precision ob- 
tained. The characteristic bump seen in the MBP and 
magnetic intergranular lane profiles is due to the Zeeman 
splitting from the high magnetic fields present in these 
components. 

2.2. Removing Oscillations 

Throughout the time-sequence, all profiles experience 
a RV shift. An examination of Fourier power in the 
pressure domain of the simulations, reveals two peaks 
at around 2 and 7 mHz. We attribute these varia- 
tions to the solar p-modes that contribute to the well- 
known '5 minute' oscillations. In order to disentangle 
the granulation signal, this oscillation needs to be re- 
moved from the profiles. Removing the oscillation allows 
the convective motions to be viewed in detail and pre- 
vents time-average component profiles from being skewed 
by oscillation-induced broadening, thus allowing a more 
accurate granulation parameterisation. 

In order to remove the effects of the oscillation, we 
shifted each component profile to fixed reference points. 
The reference points were determined by first creating 
initial time-average component profiles. The mean values 
of the line bisector wavelengths from these initial time- 
average profiles served as the reference points because 
they had no knowledge of the oscillation. We then shifted 
the components by their bisector mean wavelengths to 
preserve shape variation due to the granulation, while si- 
multaneously removing the oscillation effects. Once the 
profiles were shifted to remove the oscillation, we then 
used a final time-average for each component to be repre- 
sentative of each granulation feature. As seen in Figure^ 
there is little variation within each component category 
throughout the time-sequence, providing confidence that 
the oscillation signal does not significantly affect com- 



ponent profile shape and that it has been successfully 
removed. 
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Fig. 2. — Line profiles from the entire time series for the different 
physical components of granulation: granules (green), MBPs (red), 
magnetic (blue) and non-magnetic (purple) intergranular lanes. 
The time-average line profiles are shown in black. All components 
have been shifted by their bisector mean wavelength to the bisec- 
tor mean wavelength of their time-average profile to remove the 
oscillation signal. The breadth about each component's respective 
means is representative of the precision obtained. 



3. PARAMETERISED RECONSTRUCTION 

We examined the accuracy of the four time-average line 
profiles, representing the different components of photo- 
spheric magneto-convection, by reconstructing the line 
profiles produced by the simulations. For each granu- 
lation snapshot, we recorded the proportion of the area 
occupied by each component to the total area (filling fac- 
tor) based on the number of Stokes / profiles in a given 
category to total number in the snapshot; mean values 
for these within la are shown in Table [TJ To reconstruct 
the line profile for an entire snapshot, we multiplied each 
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Fig. 3. — The reconstructed line profiles (red) are compared to the original line profiles (black) from the 200 G simulations. Left: worst 
reconstruction achieved; right: best reconstruction achieved. 



TABLE 1 

Mean filling factor for each component with ±1ct 



Component 


Mean Filling Factor ±l<r 


Granule 


0.559 ± 0.016 


Non-Magnetic Intergranular Lane 


0.348 ± 0.019 


MBP 


0.054 ± 0.004 


Magnetic Intergranular Lane 


0.039 ± 0.006 



component profile by its corresponding filling factor and 
then added all four component profiles together. Exam- 
ples of the reconstructions are shown in Figure [3j 

We calculated the average relative error, (AI/(I)) 
(where I is the continuum intensity) , between the recon- 
structed profiles and the original profiles produced by the 
simulations. Even at the very worst, the reconstructed 
and original profiles still agree to at least ~ 2%, as seen 
in Figure [31 On average the reconstructed and original 
profiles are in agreement to ~ 0.6% or better. 

3.1. Impact of Additional Components 

We have also investigated the impact of assigning addi- 
tional categories to our parameterisation scheme. Upon 
examining the probability density distribution of the 
modulus of the vertical component of magnetic field 
(measured roughly at the continuum formation layer 
in the simulation box) versus the continuum intensity, 
shown in Figure 21 we realised a potential for an addi- 
tional two components. There are two higher density 
clusters, one at high and one at low magnetic fields with 
a range of continuum intensities. These clusters indicate 
the presence of magnetic and non-magnetic features that 
are both bright and dark. In our four component model 



we separate these regions based on a magnetic field limit 
of 1000 G and a continuum intensity limit of 0.9. How- 
ever, it is clear from Figure U that many profiles with 
intermediate magnetic fields were being included in the 
non-magnetic components (granules and non-magnetic 
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Fig. 4. — The probability density distribution of the modulus of 
the vertical component of the magnetic field approximately at the 
continuum formation level versus continuum intensity for an arbi- 
trarily chosen snapshot from the 200 G simulation. The average 
dependence of the continuum intensity on the modulus of the verti- 
cal magnetic field strength is shown by the black dot-dashed curve. 
The cutoff limits for each of the physical granulation components 
are shown in black (~ 0.9 continuum intensity, 1000 G for both 
four and six component models and, additionally, 200 G for the six 
component model). 
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Fig. 5. — Line profiles from the entire time series for the different 
physical components of granulation, as in Figure [2] with two addi- 
tional components: granules (in green), MBPs (in red), magnetic 
(in blue) and non-magnetic (in purple) intergranular lanes, semi- 
magnetic bright (in turquoise) and semi-magnetic dark (in pink) 
regions. As before, oscillation effects have been removed and the 
breadth about each component's respective means (shown in black) 
is representative of the precision obtained. 

intergranular lanes). We chose to test a model where 
there was an additional magnetic field limit at 200 G to 
separate the non-magnetic features from those with inter- 
mediate magnetic field. This additional limit then intro- 
duced two new components that we term semi-magnetic 
bright and semi-magnetic dark regions, shown in Fig- 
ure [3 with magnetic field between 200 G and 1000 G. 
We found that the six component model did not produce 
significantly smaller average relative errors when com- 
pared to the original line profiles from the simulations. 
In fact, at times the four component model even faired 
better, as seen in the top panel of Figure |6] Even at the 
most extreme, the average relative errors between four 
and six component model differed by only ~ 0.0005 (< 
1%). For simplicity we proceed with the four component 
model. 

3.2. Impact of Errors 

Ultimately we want to measure the RV signal produced 
by the convective motions integrated across a stellar disc. 
In order to do so, we must be confident in our ability to 
accurately reproduce the RVs from these convective mo- 
tions using our four component parameterisation. We 
need to ensure that any resultant RV variations are not 
due to errors introduced in the reconstruction process. 
To do so we began by measuring the RVs of the recon- 
structed profiles. 

We measured the RVs from reconstructed profiles by 
cross correlating a single, randomly selected profile with 
the entire time-series. At times, there is a possible cor- 
relation between average relative error (top panel, Fig- 
ure [6]) and RV (middle panel, Figure EJ). We believe this 
is due to the relationship between filling factor and RV. 
For example, if the granule profile has a larger error, com- 
pared to the profiles from the other components, then a 
reconstruction with a higher granule filling factor would 
have a more negative RV and also a larger error. We 
see this in Figure [6j where changes in component filling 
factor (bottom panel, shown in green for granules and 



in purple for non-magnetic intergranular lanes) correlate 
with RV (middle panel) and error (top panel). However, 
such a correlation could also indicate that the measured 
RVs may originate from errors within the reconstruction 
process rather than the granulation signal. 
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Fig. 6. — Top panel: the average relative error between the orig- 
inal line profiles from the simulations and the reconstructed line 
profiles over time. Middle panel: the RVs of the reconstructed line 
profiles over time. Four component model is shown in black and six 
component model is shown in red. Bottom panel: filling factor vs 
time for the two largest components, granules (in green) and non- 
magnetic intergranular lanes (in purple) in four-component model. 
Any correlation between RV and average relative error is likely due 
to changes in filling factors that simultaneously impact the RV and 
line shape of a profile. 

To confirm that the measured RVs for the recon- 
structed profiles do originate from the granulation sig- 
nal, and not from reconstruction errors, we compared the 
reconstructed RVs to those of the original line profiles. 
To measure the RVs of the original profiles we cross- 
correlated them with a randomly selected reconstructed 
profile. We used a reconstructed profile, arbitrarily cho- 
sen, as a reference because it has no knowledge of the 
oscillation (first panel, Figure [7} . As mentioned in Sec- 
tion l2.21 we removed the oscillation for our parameterisa- 
tion and as a result we must also remove it from the orig- 
inal profiles to directly compare with the reconstructed 
profiles. To do so, we began by estimating the oscilla- 
tion signal. If the reconstruction profiles match perfectly 
the simulated granulation then a cross-correlation of the 
original profiles with their respective reconstructed pro- 
files would provide the oscillation signal alone. While 
the reconstruction process is not perfect, we assume it is 
good enough to give a good approximation of the oscil- 
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Fig. 7. — First panel: RVs of the original profiles over the time series, the oscillation is clearly evident in the first half of the series. 
Second panel: RVs of the approximate oscillation signal. Black and red curves represent use of the four and six component reconstructions, 
respectively. Third panel: original profile RVs after approximate oscillation signal has been subtracted. Fourth panel: RVs of the 
reconstructed profiles. Fifth panel: difference between original and reconstructed RVs, after removal of oscillation. As shown in the fifth 



panel, our reconstruction is good to 
eye). 



± 20 cm s 1 on a ~100 m s 1 granulation signal (green dashed lines are plotted only to guide the 



lation signal (second panel, Figure [7J . This approximate 
oscillation signal is then subtracted from the original pro- 
file RVs, and therefore any remaining RV signal should 
be due to the convective motions and any reconstruction 
errors (third panel, Figure [7J We note here that if the os- 
cillation signal were to significantly alter the component 
profile shape, then remnant signatures would be present 
after this stage and in the reconstructions, which we do 
not find.). We then compared this remaining RV signal 
to the reconstructed RVs (fourth panel, Figure [7J , and 
calculated the difference between them to produce the 
residuals (fifth panel, Figure [JJ . The residuals are ap- 
proximately ± 20 cm s _1 for more than 95% of the time- 
sequence. As such, we conclude that our categorisation 
is good to ± 20 cm s _1 ona~ 100 m s -1 granulation 
RV variation calculated over a 12 x 12 Mm 2 area. We 
therefore conclude that any errors in the parameterisa- 
tion do not significantly impact our ability to generate 
line profiles that represent realistic granulation patterns. 

4. 50 G RECONSTRUCTION 

Our four component model parameterises granulation 
line profiles from the 200 G simulations, where the 
full spectrum of photospheric magneto-convection fea- 
tures is present. In accordance with our premise (Sec- 
tion EH]), such a parameterisation should be capable of 
accurately reproducing granulation line profiles through- 



out the quiet Sun, where average surface magnetic fields 
vary from near up to ~ 200 G. We confirm this by us- 
ing our four component 200 G parameterisation to recon- 
struct granulation profiles produced from the simulation 
with a 50 G average field. 

The 50 G simulation was parameterised in the same 
manner as the 200 G (see Section l2~Tj) . From the simula- 
tion and parameterisation, we found the filling factors of 
each of the four granulation components in a 50 G model. 
As such, we use these filling factors in conjunction with 
the four component profiles from the 200 G parameteri- 
sation to reconstruct the profiles from the 50 G simula- 
tions. As shown in Figure [51 we have achieved excellent 
accuracy in our reconstruction, ranging from at worst ~ 
1% to at best ~ 0.3%. From these results, we are con- 
fident that our parameterisation of Stokes I profiles is 
robust for the various magnetic field strengths found in 
the quiet Sun. 

5. GENERATING NEW GRANULATION LINE PROFILES 

To create new stellar line profiles that represent re- 
alistic granulation patterns over a 12 x 12 Mm 2 patch, 
we generate a probability distribution for the granules 
based on filling factors throughout the time-sequence. 
This allows us to randomly and realistically select the 
granule filling factor for a given profile. As seen in Fig- 
ure [9l there is a strong linear anti-correlation (with slope 
-1.19) between granule and non-magnetic intergranular 
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Fig. 8. — The four component line profiles from 200 G magnetic field simulation are used in conjunction with filling factors from 50 G 
simulation to reconstruct granulation line profiles in a 50 G field. The reconstructed line profiles (in red) are compared to the original line 
profiles (in black). Left: worst reconstruction achieved; right: best reconstruction achieved. 



lane filling factor (with a Pearson's correlation coefficient 
of -0.98). As a result, the granule filling factor is used in 
conjunction with this anti-correlation to determine the 
filling factor for the non-magnetic intergranular lanes. 
Together, granule and non-magnetic intergranular lanes 
make up approximately 90% of the total granulation pro- 
file (Table[T]) . Since there are no other strong correlations 
amongst filling factors, a probability distribution is again 
generated, this time for the magnetic intergranular lane 
filling factor. As the four component filling factors must 
total to one, the MBP filling factor is then determined 
by the difference between the other three filling factors 
and one. Once all four filling factors are determined we 
use these in conjunction with our four component pro- 
files determined above to create new realistic granulation 
line profiles. 

6. CONCLUDING REMARKS 

Simulating realistic solar granulation can take up to 
months for a sequence of snapshots that last little more 
than an hour. Computing the Stokes profiles from these 
simulations for just one snapshot can take hours on its 
own right. Currently, it is not feasible to tile a stellar 
disc with each 12 x 12 Mm 2 tile having an independent 
granulation pattern produced by simulations. However, 
using our parameterised model described here we can 
produce accurate granulation line profiles for magnetic 
field scenarios between 50 - 200 G. In addition, this 
method reproduces radial velocities which are accurate 
to ~ ± 20 cm s _1 on a ^100 m s _1 granulation signal. 

Throughout this paper we outlined the methodology 
for our parameterisation of granulation at disc centre. 
Naturally, the granulation component characterisation 
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Fig. 9. — There is a strong anti-correlation between granule fill- 
ing factor and the non-magnetic intergranular lane filling factor. 
The solid red curve is a linear fit. This fit is used in conjunction 
with a randomly selected granule filling factor to determine the 
non-magnetic intergranular lane filling factor used to generate new 
stellar line profiles representative of realistic granulation patterns. 

may change as we apply our parameterisation to var- 
ious disc positions, however, our methodology will re- 
main the same. Applying our parameterised method to 
off-centre regions of the stellar disc is more difficult. Pho- 
tospheric magneto-convection results in a corrugated sur- 
face. An apt analogy is to think of the g ranules as 'hills ' 
and the intergranular lanes as 'valleys' ()Dravin s 2008). 
Beyond zero degrees, we begin to see different aspects of 
these 'hills' and 'valleys.' As we incline the simulations 
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the granular walls become visible, while the intergran- 
ular lanes become ever more hidden, and MBPs disap- 
pear from view toward the limb. Consequently, there are 
changes in profile shape for granulation patterns at var- 
ious inclinations. Furthermore, higher inclinations lead 
to a physically longer line-of-sight and require finer res- 
olution (in terms of grid cell number) along the inclined 
line-of-sight due to appearance of strong under-resolved 
gradients in plasma parameters. Thus, in order to main- 
tain the same precision, an even greater amount of com- 
putational time is required. In a future paper, we will 
extend our parameterisation across the stellar disc using 
the methodology outlined here. 

Once we have a model that applies to various stel- 
lar disc positions, we will tile a stellar disc with line 
profiles that represent random, independent granulation 
patterns and integrate the wavelength-dependent inten- 
sity over the disc. These 'simulated star' observations 
will allow us to search for a way to characterise the gran- 
ulation signal as a whole by examining relationships be- 
tween line asymmetries, bisector changes and RV. As this 
method stems from the radiative 3D MHD simulations, 
we plan to expand these simulations to other spectral 
types thereby allowing our parameterised method to be 
tested beyond the solar spectrum. In future work, we will 
also explore such a parameterisation and search for noise 
removal with other stellar absorption lines-we begin here 
with one line to validate our method. There are a number 
of potential routes for achieving RV detection of Earth- 
like planets. Currently, RV measurements of exoplanets 
are done via cross-correlation with a spectral type tem- 
plate mask, requiring several thousand lines. However, 
it may be that for Earth-like confirmation we may have 
to target specific, well-behaved unblended lines, in which 
case we may only need a handle of lines modelled. Addi- 
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tionally, if we do find correlations between granular noise 
and observables then we may be able to search for said 
correlations empirically for a variety of lines, which may 
limit the need for such detailed modelling of other stellar 
lines. 

With our current parameterised four component 
model, we are able to accurately match the synthesised 
Stokes / profiles obtained from radiative diagnostics of 
3D MHD radiative simulations of the solar photosphere 
and mimic observations of solar granulation, while main- 
taining the opportunity to apply our model to other spec- 
tral types. As we move ever closer towards the confirma- 
tion of Earth-like worlds, granulation and other low-level 
noise sources s uch as meridional flows or variable gravita- 
tiona l redshift (|Beckersll2007t lMakarovll2010t iCegla et al.l 
120121 ) become an imminent limitation. While, our pri- 
mary goal is to remove/reduce granulation effects in the 
exoplanet searches, a characterisation of granulation is 
multidisciplinary and may be beneficial for solar physics 
and stellar astrophysics. Furthermore, as we expand our 
approach to other spectral types we will not only be able 
to assist exoplanet surveys around a larger variety of host 
stars, but also understand the convection process more 
thoroughly across the HR diagram. 
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